Comprehensive identification of pathogenic microbes and antimicrobial resistance genes in food products using nanopore sequencing-based metagenomics.

Foodborne pathogens, particularly antimicrobial-resistant (AMR) bacteria, remain a significant threat to global health. Given the limitations of conventional culture-based approaches, which are limited in scope and time-consuming, metagenomic sequencing of food products emerges as a promising solution. This method provides a fast and comprehensive way to detect the presence of pathogenic microbes and antimicrobial resistance genes (ARGs). Notably, nanopore long-read sequencing provides more accurate bacterial taxonomic classification in comparison to short-read sequencing. Here, we revealed the impact of food types and attributes (origin, retail place, and food processing methods) on microbial communities and the AMR profile using nanopore metagenomic sequencing. We analyzed a total of 260 food products, including raw meat, sashimi, and ready-to-eat (RTE) vegetables. Clostridium botulinum, Acinetobacter baumannii, and Vibrio parahaemolyticus were identified as the top three foodborne pathogens in raw meat and sashimi. Importantly, even with low pathogen abundance, higher percentages of samples containing carbapenem and cephalosporin resistance genes were identified in chicken and RTE vegetables, respectively. In parallel, our results demonstrated that fresh, peeled, and minced foods exhibited higher levels of pathogenic bacteria. In conclusion, this comprehensive study offers invaluable data that can contribute to food safety assessments and serve as a basis for quality indicators.

Genetic mechanisms of co-emergence of INH-resistant Mycobacterium tuberculosis strains during the standard course of antituberculosis therapy.

The incidence of isoniazid (INH) resistant Mycobacterium tuberculosis is increasing globally. This study aimed to identify the molecular mechanisms behind the development of INH resistance in M. tuberculosis strains collected from the same patients during the standard course of treatment. Three M. tuberculosis strains were collected from a patient before and during antituberculosis (anti-TB) therapy. The strains were characterized using phenotypic drug susceptibility tests, Mycobacterial Interspersed Repeated Unit-Variable-Number Tandem Repeats (MIRU-VNTR), and whole-genome sequencing (WGS) to identify mutations associated with INH resistance. To validate the role of the novel mutations in INH resistance, the mutated katG genes were electroporated into a KatG-deleted M. tuberculosis strain (GA03). Three-dimensional structures of mutated KatG were modeled to predict their impact on INH binding. The pre-treatment strain was susceptible to INH. However, two INH-resistant strains were isolated from the patient after anti-TB therapy. MIRU-VNTR and WGS revealed that the three strains were clonally identical. A missense mutation (P232L) and a nonsense mutation (Q461Stop) were identified in the katG of the two post-treatment strains, respectively. Transformation experiments showed that katG of the pre-treatment strain restored INH susceptibility in GA03, whereas the mutated katG genes from the post-treatment strains rendered negative catalase activity and INH resistance. The protein model indicated that P232L reduced INH-KatG binding affinity while Q461Stop truncated gene transcription. Our results showed that the two katG mutations, P232L and Q461Stop, accounted for the co-emergence of INH-resistant clones during anti-TB therapy. The inclusion of these mutations in the design of molecular assays could increase the diagnostic performance.IMPORTANCEThe evolution of drug-resistant strains of Mycobacterium tuberculosis within the lung lesions of a patient has a detrimental impact on treatment outcomes. This is particularly concerning for isoniazid (INH), which is the most potent first-line antimycobacterial drug. However, the precise genetic factors responsible for drug resistance in patients have not been fully elucidated, with approximately 15% of INH-resistant strains harboring unknown genetic factors. This raises concerns about the emergence of drug-resistant clones within patients, further contributing to the global epidemic of resistance. In this study, we revealed the presence of two novel katG mutations, which emerged independently due to the stress exerted by antituberculosis (anti-TB) treatment on a parental strain. Importantly, we experimentally demonstrated the functional significance of both mutations in conferring resistance to INH. Overall, this research sheds light on the genetic mechanisms underlying the evolution of INH resistance within patients and provides valuable insights for improving diagnostic performance by targeting specific mutations.

Long-Read Sequencing with Hierarchical Clustering for Antiretroviral Resistance Profiling of Mixed Human Immunodeficiency Virus Quasispecies.

BACKGROUND: HIV infections often develop drug resistance mutations (DRMs), which can increase the risk of virological failure. However, it has been difficult to determine if minor mutations occur in the same genome or in different virions using Sanger sequencing and short-read sequencing methods. Oxford Nanopore Technologies (ONT) sequencing may improve antiretroviral resistance profiling by allowing for long-read clustering. METHODS: A new ONT sequencing-based method for profiling DRMs in HIV quasispecies was developed and validated. The method used hierarchical clustering of long amplicons that cover regions associated with different types of antiretroviral drugs. A gradient series of an HIV plasmid and 2 plasma samples was prepared to validate the clustering performance. The ONT results were compared to those obtained with Sanger sequencing and Illumina sequencing in 77 HIV-positive plasma samples to evaluate the diagnostic performance. RESULTS: In the validation study, the abundance of detected quasispecies was concordant with the predicted result with the R2 of > 0.99. During the diagnostic evaluation, 59/77 samples were successfully sequenced for DRMs. Among 18 failed samples, 17 were below the limit of detection of 303.9 copies/µL. Based on the receiver operating characteristic analysis, the ONT workflow achieved an F1 score of 0.96 with a cutoff of 0.4 variant allele frequency. Four cases were found to have quasispecies with DRMs, in which 2 harbored quasispecies with more than one class of DRMs. Treatment modifications were recommended for these cases. CONCLUSIONS: Long-read sequencing coupled with hierarchical clustering could differentiate the quasispecies resistance profiles in HIV-infected samples, providing a clearer picture for medical care.

Identification of microbial community in the urban environment: The concordance between conventional culture and nanopore 16S rRNA sequencing.

Introduction: Microbes in the built environment have been implicated as a source of infectious diseases. Bacterial culture is the standard method for assessing the risk of exposure to pathogens in urban environments, but this method only accounts for <1% of the diversity of bacteria. Recently, full-length 16S rRNA gene analysis using nanopore sequencing has been applied for microbial evaluations, resulting in a rise in the development of long-read taxonomic tools for species-level classification. Regarding their comparative performance, there is, however, a lack of information. Methods: Here, we aim to analyze the concordance of the microbial community in the urban environment inferred by multiple taxonomic classifiers, including ARGpore2, Emu, Kraken2/Bracken and NanoCLUST, using our 16S-nanopore dataset generated by MegaBLAST, as well as assess their abilities to identify culturable species based on the conventional culture results. Results: According to our results, NanoCLUST was preferred for 16S microbial profiling because it had a high concordance of dominant species and a similar microbial profile to MegaBLAST, whereas Kraken2/Bracken, which had similar clustering results as NanoCLUST, was also desirable. Second, for culturable species identification, Emu with the highest accuracy (81.2%) and F1 score (29%) for the detection of culturable species was suggested. Discussion: In addition to generating datasets in complex communities for future benchmarking studies, our comprehensive evaluation of the taxonomic classifiers offers recommendations for ongoing microbial community research, particularly for complex communities using nanopore 16S rRNA sequencing.

Probable airborne transmission of Burkholderia pseudomallei causing an urban outbreak of melioidosis during typhoon season in Hong Kong, China.

Between January 2015 and October 2022, 38 patients with culture-confirmed melioidosis were identified in the Kowloon West (KW) Region, Hong Kong. Notably, 30 of them were clustered in the Sham Shui Po (SSP) district, which covers an estimated area of 2.5 km2. Between August and October 2022, 18 patients were identified in this district after heavy rainfall and typhoons. The sudden upsurge in cases prompted an environmental investigation, which involved collecting 20 air samples and 72 soil samples from residential areas near the patients. A viable isolate of Burkholderia pseudomallei was obtained from an air sample collected at a building site five days after a typhoon. B. pseudomallei DNA was also detected in 21 soil samples collected from the building site and adjacent gardening areas using full-length 16S rRNA gene sequencing, suggesting that B. psuedomallei is widely distributed in the soil environment surrounding the district. Core genome-multilocus sequence typing showed that the air sample isolate was phylogenetically clustered with the outbreak isolates in KW Region. Multispectral satellite imagery revealed a continuous reduction in vegetation region in SSP district by 162,255 m2 from 2016 to 2022, supporting the hypothesis of inhalation of aerosols from the contaminated soil as the transmission route of melioidosis during extreme weather events. This is because the bacteria in unvegetated soil are more easily spread by winds. In consistent with inhalational melioidosis, 24 (63.2%) patients had pneumonia. Clinicians should be aware of melioidosis during typhoon season and initiate appropriate investigation and treatment for patients with compatible symptoms.

Evaluation of Mycobacterium tuberculosis enrichment in metagenomic samples using ONT adaptive sequencing and amplicon sequencing for identification and variant calling.

Sensitive detection of Mycobacterium tuberculosis (TB) in small percentages in metagenomic samples is essential for microbial classification and drug resistance prediction. However, traditional methods, such as bacterial culture and microscopy, are time-consuming and sometimes have limited TB detection sensitivity. Oxford nanopore technologies (ONT) MinION sequencing allows rapid and simple sample preparation for sequencing. Its recently developed adaptive sequencing selects reads from targets while allowing real-time base-calling to achieve sequence enrichment or depletion during sequencing. Another common enrichment method is PCR amplification of the target TB genes. In this study, we compared both methods using ONT MinION sequencing for TB detection and variant calling in metagenomic samples using both simulation runs and those with synthetic and patient samples. We found that both methods effectively enrich TB reads from a high percentage of human (95%) and other microbial DNA. Adaptive sequencing with readfish and UNCALLDE achieved a 3.9-fold and 2.2-fold enrichment compared to the control run. We provide a simple automatic analysis framework to support the detection of TB for clinical use, openly available at https://github.com/HKU-BAL/ONT-TB-NF . Depending on the patient's medical condition and sample type, we recommend users evaluate and optimize their workflow for different clinical specimens to improve the detection limit.

The clinical utility of Nanopore 16S rRNA gene sequencing for direct bacterial identification in normally sterile body fluids.

The prolonged incubation period of traditional culture methods leads to a delay in diagnosing invasive infections. Nanopore 16S rRNA gene sequencing (Nanopore 16S) offers a potential rapid diagnostic approach for directly identifying bacteria in infected body fluids. To evaluate the clinical utility of Nanopore 16S, we conducted a study involving the collection and sequencing of 128 monomicrobial samples, 65 polymicrobial samples, and 20 culture-negative body fluids. To minimize classification bias, taxonomic classification was performed using 3 analysis pipelines: Epi2me, Emu, and NanoCLUST. The result was compared to the culture references. The limit of detection of Nanopore 16S was also determined using simulated bacteremic blood samples. Among the three classifiers, Emu demonstrated the highest concordance with the culture results. It correctly identified the taxon of 125 (97.7%) of the 128 monomicrobial samples, compared to 109 (85.2%) for Epi2me and 102 (79.7%) for NanoCLUST. For the 230 cultured species in the 65 polymicrobial samples, Emu correctly identified 188 (81.7%) cultured species, compared to 174 (75.7%) for Epi2me and 125 (54.3%) for NanoCLUST. Through ROC analysis on the monomicrobial samples, we determined a threshold of relative abundance at 0.058 for distinguishing potential pathogens from background in Nanopore 16S. Applying this threshold resulted in the identification of 107 (83.6%), 117 (91.4%), and 114 (91.2%) correctly detected samples for Epi2me, Emu, and NanoCLUST, respectively, in the monomicrobial samples. Nanopore 16S coupled with Epi2me could provide preliminary results within 6 h. However, the ROC analysis of polymicrobial samples exhibited a random-like performance, making it difficult to establish a threshold. The overall limit of detection for Nanopore 16S was found to be about 90 CFU/ml.

Clinical utility of target amplicon sequencing test for rapid diagnosis of drug-resistant Mycobacterium tuberculosis from respiratory specimens.

An in-house-developed target amplicon sequencing by next-generation sequencing technology (TB-NGS) enables simultaneous detection of resistance-related mutations in Mycobacterium tuberculosis (MTB) against 8 anti-tuberculosis drug classes. In this multi-center study, we investigated the clinical utility of incorporating TB-NGS for rapid drug-resistant MTB detection in high endemic regions in southeast China. From January 2018 to November 2019, 4,047 respiratory specimens were available from patients suffering lower respiratory tract infections in Hong Kong and Guangzhou, among which 501 were TB-positive as detected by in-house IS6110-qPCR assay with diagnostic sensitivity and specificity of 97.9 and 99.2%, respectively. Preliminary resistance screening by GenoType MTBDRplus and MTBDRsl identified 25 drug-resistant specimens including 10 multidrug-resistant TB. TB-NGS was performed using MiSeq on all drug-resistant specimens alongside 67 pan-susceptible specimens, and demonstrated 100% concordance to phenotypic drug susceptibility test. All phenotypically resistant specimens with dominating resistance-related mutations exhibited a mutation frequency of over 60%. Three quasispecies were identified with mutation frequency of less than 35% among phenotypically susceptible specimens. They were well distinguished from phenotypically resistant cases and thus would not complicate TB-NGS results interpretations. This is the first large-scale study that explored the use of laboratory-developed NGS platforms for rapid TB diagnosis. By incorporating TB-NGS with our proposed diagnostic algorithm, the workflow would provide a user-friendly, cost-effective routine diagnostic solution for complicated TB cases with an average turnaround time of 6 working days. This is critical for timely management of drug resistant TB patients and expediting public health control on the emergence of drug-resistant TB.

A multi-omics investigation into the mechanisms of hyper-virulence in Mycobacterium tuberculosis.

Clinical manifestations of tuberculosis range from asymptomatic infection to a life-threatening disease such as tuberculous meningitis (TBM). Recent studies showed that the spectrum of disease severity could be related to genetic diversity among clinical strains of Mycobacterium tuberculosis (Mtb). Certain strains are reported to preferentially invade the central nervous system, thus earning the label "hypervirulent strains".However, specific genetic mutations that accounted for enhanced mycobacterial virulence are still unknown. We previously identified a set of 17 mutations in a hypervirulent Mtb strain that was from TBM patient and exhibited significantly better intracellular survivability. These mutations were also commonly shared by a cluster of globally circulating hyper-virulent strains. Here, we aimed to validate the impact of these hypervirulent-specific mutations on the dysregulation of gene networks associated with virulence in Mtb via multi-omic analysis. We surveyed transcriptomic and proteomic differences between the hyper-virulent and low-virulent strains using RNA-sequencing and label-free quantitative LC-MS/MS approach, respectively. We identified 25 genes consistently differentially expressed between the strains at both transcript and protein level, regardless the strains were growing in a nutrient-rich or a physiologically relevant multi-stress condition (acidic pH, limited nutrients, nitrosative stress, and hypoxia). Based on integrated genomic-transcriptomic and proteomic comparisons, the hypervirulent-specific mutations in FadE5 (g. 295,746 C >T), Rv0178 (p. asp150glu), higB (p. asp30glu), and pip (IS6110-insertion) were linked to deregulated expression of the respective genes and their functionally downstream regulons. The result validated the connections between mutations, gene expression, and mycobacterial pathogenicity, and identified new possible virulence-associated pathways in Mtb.

Screening Repurposed Antiviral Small Molecules as Antimycobacterial Compounds by a Lux-Based phoP Promoter-Reporter Platform.

The emergence of multidrug-resistant strains and hyper-virulent strains of Mycobacterium tuberculosis are big therapeutic challenges for tuberculosis (TB) control. Repurposing bioactive small-molecule compounds has recently become a new therapeutic approach against TB. This study aimed to identify novel anti-TB agents from a library of small-molecule compounds via a rapid screening system. A total of 320 small-molecule compounds were used to screen for their ability to suppress the expression of a key virulence gene, phop, of the M. tuberculosis complex using luminescence (lux)-based promoter-reporter platforms. The minimum inhibitory and bactericidal concentrations on drug-resistant M. tuberculosis and cytotoxicity to human macrophages were determined. RNA sequencing (RNA-seq) was conducted to determine the drug mechanisms of the selected compounds as novel antibiotics or anti-virulent agents against the M. tuberculosis complex. The results showed that six compounds displayed bactericidal activity against M. bovis BCG, of which Ebselen demonstrated the lowest cytotoxicity to macrophages and was considered as a potential antibiotic for TB. Another ten compounds did not inhibit the in vitro growth of the M. tuberculosis complex and six of them downregulated the expression of phoP/R significantly. Of these, ST-193 and ST-193 (hydrochloride) showed low cytotoxicity and were suggested to be potential anti-virulence agents for M. tuberculosis.

The Clinical Utility of Two High-Throughput 16S rRNA Gene Sequencing Workflows for Taxonomic Assignment of Unidentifiable Bacterial Pathogens in Matrix-Assisted Laser Desorption Ionization-Time of Flight Mass Spectrometry.

Bacterial pathogens that cannot be identified using matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) are occasionally encountered in clinical laboratories. The 16S rRNA gene is often used for sequence-based analysis to identify these bacterial species. Nevertheless, traditional Sanger sequencing is laborious, time-consuming, and low throughput. Here, we compared two commercially available 16S rRNA gene sequencing tests that are based on Illumina and Nanopore sequencing technologies, respectively, in their ability to identify the species of 172 clinical isolates that failed to be identified by MALDI-TOF MS. Sequencing data were analyzed by the respective built-in programs (MiSeq Reporter software of Illumina and Epi2me of Nanopore) and BLAST+ (v2.11.0). Their agreement with Sanger sequencing on species-level identification was determined. Discrepancies were resolved by whole-genome sequencing. The diagnostic accuracy of each workflow was determined using the composite sequencing result as the reference standard. Despite the high base-calling accuracy of Illumina sequencing, we demonstrated that the Nanopore workflow had a higher taxonomic resolution at the species level. Using built-in analysis algorithms, the concordance of Sanger 16S with the Illumina and Nanopore workflows was 33.14% and 87.79%, respectively. The agreement was 65.70% and 83.14%, respectively, when BLAST+ was used for analysis. Compared with the reference standard, the diagnostic accuracy of Nanopore 16S was 96.36%, which was identical to that of Sanger 16S and better than that of Illumina 16S (69.07%). The turnaround time of the Illumina workflow and the Nanopore workflow was 78 h and 8.25 h, respectively. The per-sample cost of the Illumina and Nanopore workflows was US$28.5 and US$17.7, respectively.

A low-cost TaqMan minor groove binder probe-based one-step RT-qPCR assay for rapid identification of N501Y variants of SARS-CoV-2.

The increasing prevalence of N501Y variants of SARS-CoV-2 has kindled global concern due to their enhanced transmissibility. Genome sequencing is the gold standard method to identify the emerging variants of concern. But it is time-consuming and expensive, limiting the widespread deployment of genome surveillance in some countries. Health authorities surge the development of alternative assay to expand screening capacity with reduced time and cost. In this study, we developed an in-house TaqMan minor groove binder (MGB) probe-based one-step RT-qPCR assay to detect the presence of N501Y mutation in SARS-CoV-2. A total of 168 SARS-CoV-2 positive respiratory specimens were collected to determine diagnostic accuracy of the RT-qPCR assay. As a reference standard, PANGO lineages and the mutation patterns of all samples were characterised by whole-genome sequencing. The analytical sensitivity and the ability of the assay to detect low frequency of N501Y variants were also evaluated. A total of 31 PANGO lineages were identified from 168 SARS-CoV-2 positive cases, in which 34 samples belonged to N501Y variants, including B.1.1.7 (n = 20), B.1.351 (n = 12) and P.3 (n = 2). The N501Y RT-qPCR correctly identified all 34 samples as N501Y-positive and the other 134 samples as wildtype. The limit-of-detection of the assay consistently achieved 1.5 copies/µL on four different qPCR platforms. N501Y mutation was successfully detected at an allele frequency as low as 10 % in a sample with mixed SARS-CoV-2 lineage. The N501Y RT-qPCR is simple and inexpensive (US$1.6 per sample). It enables robust high-throughput screening for surveillance of SARS-CoV-2 variants of concern harbouring N501Y mutation.

Territorywide Study of Early Coronavirus Disease Outbreak, Hong Kong, China.

Initial cases of coronavirus disease in Hong Kong were imported from mainland China. A dramatic increase in case numbers was seen in February 2020. Most case-patients had no recent travel history, suggesting the presence of transmission chains in the local community. We collected demographic, clinical, and epidemiologic data from 50 patients, who accounted for 53.8% of total reported case-patients as of February 28, 2020. We performed whole-genome sequencing to determine phylogenetic relationship and transmission dynamics of severe acute respiratory syndrome coronavirus 2 infections. By using phylogenetic analysis, we attributed the community outbreak to 2 lineages; 1 harbored a common mutation, Orf3a-G251V, and accounted for 88.0% of the cases in our study. The estimated time to the most recent common ancestor of local coronavirus disease outbreak was December 24, 2019, with an evolutionary rate of 3.04 × 10-3 substitutions/site/year. The reproduction number was 1.84, indicating ongoing community spread.

Targeted-Sequencing Workflows for Comprehensive Drug Resistance Profiling of Mycobacterium tuberculosis Cultures Using Two Commercial Sequencing Platforms: Comparison of Analytical and Diagnostic Performance, Turnaround Time, and Cost.

BACKGROUND: The emergence of Mycobacterium tuberculosis with complex drug resistance profiles necessitates a rapid and comprehensive drug susceptibility test for guidance of patient treatment. We developed two targeted-sequencing workflows based on Illumina MiSeq and Nanopore MinION for the prediction of drug resistance in M. tuberculosis toward 12 antibiotics. METHODS: A total of 163 M. tuberculosis isolates collected from Hong Kong and Ethiopia were subjected to a multiplex PCR for simultaneous amplification of 19 drug resistance-associated genetic regions. The amplicons were then barcoded and sequenced in parallel on MiSeq and MinION in respective batch sizes of 24 and 12 samples. A web-based bioinformatics pipeline, BacterioChek-TB, was developed to translate the raw datasets into clinician-friendly reports. RESULTS: Both platforms successfully sequenced all samples with mean read depths of 1,127× and 1,649×, respectively. The variant calling by MiSeq and MinION could achieve 100% agreement if variants with an allele frequency of <40% reported by MinION were excluded. Both workflows achieved a mean clinical sensitivity of 94.8% and clinical specificity of 98.0% when compared with phenotypic drug susceptibility test (pDST). Turnaround times for the MiSeq and MinION workflows were 38 and 15 h, facilitating the delivery of treatment guidance at least 17-18 days earlier than pDST, respectively. The higher cost per sample on the MinION platform ($71.56) versus the MiSeq platform ($67.83) was attributed to differences in batching capabilities. CONCLUSION: Our study demonstrates the interchangeability of MiSeq and MinION platforms for generation of accurate and actionable results for the treatment of tuberculosis.